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Abstract 

An engineering approach to the design of an in-line dynamic absorber is presented using 
transfer matrix techniques for an infinitely long cable system. The natural frequencies of 
the system are determined and the properties of the installed damper are optimized to 
maximize the amount of energy absorbed. Numerical simulations show that the optimally 
tuned and damped in-line dynamic absorber will dissipate 50% of incident transverse 
vibration energy. A design methodology for a nominal marine monitoring program is 
presented and specific engineering constraints provided to ensure that the modeled system 
meets design specifications 
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Chapter 1 



Introduction 

Investigation of the marine environment is a growing industry. In order to 
economically obtain enough information, increased usage of remote sensing technologies 
will be required. The use of these technologies will require that such information 
collection systems be robust enough to survive in the harsh environment for prolonged 
periods of time without human intervention. In most cases, the areas of concern for the 
scientist can be defined by a geographical boundary, and an ocean mooring can be used to 
locate the measuring devices. The motions of the moored body can damage the 
equipment or produce erroneous readings. Thus, the robust requirement of the 
measurement system will lead to either a hardening of the sensitive electronics of the 
system or other some other sort of isolation method that will decrease the cyclic loading 
that will be seen by the system. 

The problem of the motion of a moored buoy with sensitive electronic instruments 
located along the length of the cable can be modeled as a string with a rigid lumps. A 
previous analytical approach of the problem was performed by Li [1]. In his analysis, Li 
was able to provide a analytical description of how to model a string/rigid lump system by 
the geometric and material properties of the components and the cable and consider the 
system more than just a lumped mass at a particular point along the cable. The energy of 
the wave traveling along the cable might be transmitted through the component or 
reflected back, depending upon the specifics of the problem. 

If the physical connection of the equipment along the length of the cable were to 
include a spring and a damper, the model becomes more complex. The energy of the 
wave traveling along the cable might then be absorbed, transmitted or reflected at the 
interface, dependent upon the nature of the system and the frequency of the incoming 
wave. Li was able to investigate the modeling of a mass/string/damper system that would 
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be installed along the line, as well as endpoint terminations of the cable — either connected 
to a moving or a fixed frame of reference. The in-line system has more of a practical 
application as it can become a part of the device itself, as components may be changed 
along the length of a moored cable in order to alter the measurement systems available to 
the scientist in the field. 

The idea that the in-line component could be isolated via an appropriately 
engineered system has significant applications. The component then could be designed 
more towards the gathering of information, not just being able to withstand the normal 
marine forces. Using the ideas established by Li, this thesis will investigate the 
characteristics of such a moored body system and develop an engineering approach to 
determining the optimum value of the constituent parts of the measurement devices such 
that the operation of the in-line components is not impaired. 

This thesis will investigate the design of an in-line dynamic absorber tuned to the 
frequency where the maximum amount of energy traveling along the ocean mooring 
system can be absorbed. Chapter 2 outlines the rationale that drives the need for this 
analysis. Additionally, an application that presents a process for the design of a marine 
monitoring program and establishes technical requirements is introduced. 

Chapter 3 is the introduction to the use of transfer matrix notation and the benefits 
that result in the use of this system. Chapter 4 is the analysis section of this thesis. The 
search for an analytical approach to determining the optimum amount of energy that can 
be absorbed in the mooring system is provided. A numerical result is obtained and the 
optimum value for the absorber for a given system is provided. Chapter 5 is a discussion 
of the major findings and potential areas for further research. 
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Chapter 2 



Oceanographic Monitoring System 

2.1 Introduction 

The design process to implement an oceanographic monitoring system discussed in 
chapters 1 must minimize the effects of marine induced forces on an infinitely long cable 
system with installed monitoring equipment. One practical application of the design 
process will involve the use of a monitoring system that is located in a marine environment 
of interest. This chapter will describe how such a monitoring plan is designed and 
implemented. 

The first step in the development of a monitoring plan requires an executive level 
design methodology. Using the concepts from a report by the National Research Council 
[2], the following steps are used: 

1 . Define Expectations and Goals 

2. Define Study Strategy 

3. Conduct Exploratory Studies if Needed 

4. Develop Sampling Design 

5. Implement Study 

The expectations and goals of the particular monitoring plan would be based upon 
the marine site chosen. The ultimate goal is to have a monitoring system that produces 
useful information that can used to make a reasoned decision. Thus, the effects that are 
seen in the monitoring results must reflect information about the situation that exists in the 
marine environment. The legal framework that guides the administration of human actions 
in the marine environment also molds the expectations and goals of the monitoring plan. 

Defining the study strategy of the monitoring plan requires the implementation of a 
strategy that begins with the identification of resources in the marine environment that are 
at risk. From that identification, a conceptual model is developed that establishes links 
that exist in the marine ecosystem components. These links involve spatial, temporal, 
biological and chemical variables that become the basis of the study strategy. Feedback 
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among the monitored variables in the plan might dictate that some be changed to better 
reflect the marine environment. The inherent uncertainties of a conceptual model must 
also be considered in the definition of the study strategy. The final result is a set of 
questions concerning the marine environment that can be tested to describe the state of the 
marine environment. 

Exploratory studies that have been performed on similar marine ecosystems can be 
used as a basis for the monitoring plan. If the need exists to validate certain technical 
issues such as the measurement of parameters under various hydrodynamic conditions (for 
a coastal system, high river flow from spring melt of snow or low river flow from a period 
of no precipitation ), preliminary studies should be performed. The primary research 
performed for this validation becomes an interdependent part of the monitoring plan, 
altered as the need arises to obtain specific information. 

The development of a sampling plan uses the testable questions described 
previously in the study plan and identifies meaningful levels of changes that can be 
measured. In the establishment of these meaningful levels of change, the natural variability 
of the marine environment in question needs to be considered. A meaningful level of 
change must also be established that will be observable in a time frame that allows a 
reasoned decision to be made that will have an impact on the system being measured. This 
definition of meaningful change will also be as fluid as the improvement of prediction 
models and the analysis of data will allow. The selection of the variables to measure in the 
monitoring plan is limited by the resources available, and they must be chosen to provide 
the most useful information about the status of the system. Thus, variables may be 
measured that will provide a realistic view of the system, although they may not be part of 
the predictive model. These substitute variables may contain early warning signs of 
system degradation or information concerning specific items of interest to the users of the 
sampling plan. In either case, the substitute variables will provide information that can be 
applied to the predictive model via knowledge of the interactions between the substitute 
variables and the model variables. The use of specific oceanographic measurement 
equipment is then chosen to reflect the variables measured. 
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Two concerns that must be addressed in the design of the sampling plan are the 
level of confidence in the knowledge of the interactions among the elements of the system 
and the statistical characteristics of the sampling of the designated variables. If the 
interactions among the biological elements is purely theory, then the uncertainty must be 
addressed in the analysis. If the sampling process requires a significant number of samples 
in order to verify the validity of the measurement, then the monitoring plan is not efficient 
in obtaining the maximum amount of information for the effort. 

Using the concepts described in the preceding paragraphs, the monitoring plan will 
be implemented, with constant feedback to ensure that the model, the sampling plan and 
the uncertainty of the monitoring plan are all considered when updating the process. 

2.2 Tiered System Structure 

Using the guidelines for establishing a monitoring plan described above leads to a 
variety of monitoring plans based upon conditional aspects of the system being studied. 

An approach used by a team of graduate students at MIT [9]. provided a tiered system 
approach to a monitoring plan that allows a manager to best allocate the available 
resources when developing the monitoring plan. 




Figure 2- 1: Three Tiered Strategy 
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The tiered structure is classified to characterize three potential areas of impact 
within a monitoring plan area. Those areas are: 1) Water Column Region, 2) Living 
Resources , and 3) Sediment/Benthos. The monitoring plan addresses these three areas in 
three tiered levels of measurement. The three tiered strategy is graphically represented in 
figure 4-1. Each box represents a variable that is part of the generic design study. The 
number of the boxes would change based upon the specific application The arrows 
represent that the tiers are logically linked to each other. Note that the logical links from 
one tier to the next may indicate more than one cause and effect relationship for a 
particular variable. 

The first tier in the monitoring structure is comprised of variables that react 
quickly to the changes in the marine system, and provide advanced warning of potential 
events. Tier I items include the measurement of water quality variables, such as 
temperature, salinity and dissolved oxygen in the water column. The living resources 
classification is measured in fish population studies for tier I. The benthos region tier I 
measurements include bottom observations that may be made with remote sensing 
equipment. 

If the variables measured in the first tier exceed a threshold value, then more 
extensive analysis must be performed. This triggers the investigation of the tier II 
variables via a logical link to the most probable source of the penurbation. The tier II 
variables require a more complex analysis that is generally more expensive and not be 
performed in situ. An example of a tier II variable is the chemical analysis of the water 
column that requires a grab sample be obtained and transported to a laboratory. The 
living resources tier II variable would be the analysis of animal or plant material to 
determine its make-up and compare it to standard information. A tier II benthos region 
variable would again be a chemical analysis of the material of the sediment. 

The final portion of the tiered system is the logical link from the identification of 
the perturbation to the most probable source of the change. This involves the tracing of 
the cause of the perturbation to its source. The basis for this would be set up in the design 
study, when potential sources of perturbation to the marine system would be identified as 
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any facility that could potentially add material and or heat to the region around the 
boundary of the marine system being monitored (power plants, industrial plants, chemical 
plants). This detective effort is the final result of the monitoring plan because it is the 
information that managers need to form their reasoned opinions concerning liability and 
potential litigation. 

The most significant advantage of the three tiered system is its modularity. The 
system can promptly be updated to reflect the information obtained when the monitoring 
system goes into effect. Additionally, the three classifications of water column, living 
resources and benthos/sediment characterize the essential aspects of the marine 
environment. Each tier can be expanded horizontally as information is obtained and more 
classifications can be added if required. If additional links are provided with the 
improvement in knowledge, the links can also be updated. 

The modularity of the system allows the resource manager to best use the assets 
that are available if the entire plan can not be fully implemented. The logical links that 
connect the tiers describe flow paths that can be chosen in either a bottom up or top down 
approach to investigate the potential cause and effect relationships. 

2.3 Implementation 

The use of moored oceanographic systems to measure components that are part of 
the monitoring plan will likely be used because of the cost effectiveness of remote sensing 
and the possibility of real time data management from a shore based office. In order to 
technically achieve the use of moored oceanographic systems, it is imperative that external 
forcing functions that move the system be controlled to allow sensitive equipment to 
operate without damage. 

The use of the design process to design a dynamic absorber can be used most 
effectively in the tier I analysis of the monitoring plan. The vision is that the particular 
measurement desired to be measured in the water column m situ can be designed to 
operate installed on a moored buoy system. The power source for the system can be 
designed to last for a designated period of time that the user would correspond to a 
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maintenance schedule for the equipment. Current renewable energy systems rely on solar 
energy or other means that would be installed in the moored system near the free surface 
that would suffer if excess drag or vibrations caused the buoyant portion of the moored 
system to go underwater. Communications with the land based researchers would also be 
connected to the installed devices via the mooring cable. This communication suite would 
most likely be installed in a buoyant body near the free surface in a vulnerable position. 

A specific example of a typical oceanographic system would include acoustic 
transducers to monitor various components of the marine environment. These transducers 
are very sensitive to excessive vibrations, and the design process to remove as much 
vibration energy from the system as possible would help to ensure that the installed 
equipment could operate properly. 

With the implementation of Autonomous Underwater Vehicles (AUVs) in the 
marine environment, certain tier II variables can also be measured. Ongoing work within 
the Ocean Engineering Department at MIT is investigating the implementation of new 
technology that will greatly reduce the amount of power required to operate the 
propulsion systems of AUVs. The computer program required in order to operate an 
AUV in the underwater marine environment has also been the subject of recent studies at 
MIT. The power and communication source for the AUV would be connected to another 
moored cable system that would have an underwater docking port. This docking port 
would provide the communications link with the land based researchers and replenish the 
power used in the tier II investigations. 

Of particular concern to each of the applications described above is the principal of 
flow induced vibrations. Even with a steady current moving past a cable system or a cable 
system being towed by a vessel at a constant speed, oscillatory motions are induced when 
vortex shedding occurs on the cable. The motions of the cable are dominated by motions 
that are transverse to the direction of the flow and occur at a frequency equal to the vortex 
shedding frequency. In-line motions occur at twice the vortex shedding frequency, and 
are not considered in this analysis because of their small contribution to the cable motion. 
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The vortex shedding excitation frequency often coincides with a resonant frequency of the 
cable and results in substantial response. 

The non-dimensional parameter that characterizes the frequency of vortex 
shedding is the Strouhal number, and the typical value associated with it is .21 and it is 
constant over a wide range of Reynolds number. [4]. The expression relating the Strouhal 
number and the vortex shedding frequency is: 



EQN 2- 1 

Where to is the angular frequency, V is the velocity of the cable being towed or the current 
flowing past the cable, and d is the diameter of the cable. In the design process, a limit 
needs to be established for reasonable expectations of vortex shedding frequencies that the 
moored system may face. Assuming a limit of a 5 knot current (2.54 m/sec) and a range 
of cable diameters from 1/4 inch to 1 inch. Figure 2-2 represents the expected vortex 
shedding frequencies. 




Figure 2- 2: Vortex Shedding Excitation Frequencies in a 2.5 m/sec Current 

In addition to this vibrational lock in of the cable, an additional concern for flow 
induced vibration is the increase in the nondimensional drag coefficient, CD. This increase 
in drag has major implication to mooring design. 
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Understanding the need for the oceanographic monitoring and the environment in 
which the components will be operated leads directly to an investigation to reduce the 
amount of energy that the installed components will see. This investigation to limit the 
effect of the ocean loading, in particular the vortex-induced vibrations, is pursued using an 
in-line dynamic absorber. 
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Chapter 3 



Transfer Matrix Evaluation 

3.1 Background 

The transfer matrix method is a technique to evaluate a complicated system 
through the use of discrete transformations from one location to another. For a sequential 
system of cable, rigid bodies, springs and dampers, the final system transfer matrix is the 
product of each transfer matrix from one location to the other. The benefit of this system 
is the ability to evaluate the response of a complicated system through sequential 
multiplication of finite transfer matrices from point to point. The analysis of this thesis 
will evaluate the transfer matrix for an oceanographic moored system with a damping 
connector installed in-line, shown in Figure 4-1. 

Using the transfer matrix method requires that a state vector be defined. The state 
vector is comprised of Fourier transforms of displacements and the corresponding internal 
forces [5]. In the case of a cable system with installed masses, springs and dampers, the 
state vector is described by: 




EQN 3- 1 



where y is the Fourier transform of the displacement of the system, y' is the Fourier 



transform of the rotation of the system, and T is the tension of the system. .AJl three items 
are evaluated at the particular point (annotated by the subscript /), and the transformation 
from one location to the next is defined by: 



= 



EQN 3- 2 



where the transfer matrix links the two state vectors. 
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3.2 Derivation of a Transfer Matrix 



The determination of the transfer matrix is based upon an analytical description of 
the cable system. This thesis uses the transfer matrix library determined by Li.[l], In 
order to describe how the transfer matrices are established, an alternate description of an 
equivalent system is provided. Knowing the geometric relationships and determining the 
equations of motion are the keys to deriving the transfer matrix for a system. 




Figure 3-1 : Schematic of System 



Using figure 3-1 and setting Mi=M 2 =M, the following equations of motion can be 
described in the y direction; 

^yc2 = T,yi- t ; 

+A<Vc. = T,.,yU EQN 3- 3 

^CPc2 + v„) = Ty-T,_y,_^ 



Where 7^ is the internal force at the joint located at yo acting in the yo direction 

The equations for summation of moments about the center of gravity of the two 
links are described by the following; 
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EQN 3- 4 



-e,) + Ud, -9,)+ T,{y: -, 9,)f 
I A = -9,) + R,{9, - ^, ) + 7 ;., Cy;_, - 9 , )d EQN 3- 5 



Where K« is the value of the torsional spring constant (Nm/rad) and is the value of the 

torsional damping constant (Nm/(rad/sec)). 

If the following geometric relations: 



c d 

>'cl = —/,-! + —>'<, 
Ly Ly 



f ^ 



EQN 3- 6 
EQN 3- 7 



are inserted in Equation 3-3, we obtain; 



= 7tA7T7(?:X 3^.-1 



i.d+f)M 






EQN 3- 8 



And if we substitute the following angular relationships; 

and 

A ■ -4 



EQN 3- 9 



into equations 3-4 and 3-5, we obtain; 

, (ZzM = X. + Ji^^yozju . AzAy + 7 -(y . 2IZI1) j 
L, L, u L, U 4 



EQN 3- 10 



and 



j >0 - j>... ^ + 7 ;., Cy;., - y 



4 



4 4 



4 4 



4 

EQN 3- 11 
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Performing the Fourier transform of Equation 3-8 and setting T,=Ti.i=T, d=f=l 2 , c=e=li, 
and Li=L 2 =L, we obtain: 

= -- T^b'. +y.) EQN 3- 12 

0 ) J 2/, 

Where , and y\ are the Fourier transforms of the displacements and 

the slopes at the particular locations. 

Performing the Fourier transform of Equation 3-10 and Ic2 = Mk^ (k is the radius 
of gyration), we obtain: • 



. K, -O’.-. EQN 3- 13 

Z/ L Lf Z/ 

Performing the Fourier transform of Equation 3-1 1, and setting Id =’ Mk^ (k is the radius 
of gyration), we obtain: 



Lf Lf L 0 L 0 

14 



If we substitute Equation 3-12 into Equation 3-13 and rewrite, we obtain: 






1-1 



EQN 3- 15 



Where: 



, -MvV‘, /, , K, I, Kg jcoRg l.Rgjo) LT IJ 
A = (1 + _!_)+_! + _!_«+ 2 Z + _L_22_ + _=_ + _L_ 

' L IL L LI. L LL L 2L 



B, = 



-Mk^T TKg RgT 



2/2 Lq)‘ jo)l^ 



-TL + 



' 2o)' 



^ Mk CO ^ Kg I, Kg j(oR„ l,R„jco 
' 2Ll^ ' L Ll^ L LI, 
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^ _ Mk-T TKg R^T Tl, , T‘ 

' ” 24 Lq)‘ jcoU 2L ^ 2o)- 



And combining Equations 3-13 and 3-14 and rewriting to obtain: 



^2y>+B2y, +D,y,_, 



EQN 3- 16 



Where 

, -MkW ,Tl, ^ ^ ^ ^ -MkW TL ^ ^ „ 

^ L L * ' - L L ‘ ‘ 



The final equations for the transfer function can be formed using the following 
relationship: 



1 

1 

1 




'4i 


^2 


1 

V 1 
7 

1 













Thus: 






1 -B,Dy 

T A,B^^ - AyB, 



^ _ j, A,C^ - AX, 

■' ~ A,B, - A.B, 



A,D^ - A^D, 

AjB^ — A^Bf 



EQN 3- 17 



Using the methods of transfer matrix operations allow an efficient calculation of 
important parameters such as the amount of energy that is reflected, absorbed or 
transmitted through a system. The analysis in the next chapter will provide that 
information and the techniques needed in order to reduce the amount of flow-induced 
vibrations that are seen in the system 
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Chapter 4 



In-line Dynamic Absorber 

4.1 System Description 

Using the principles described in chapter 3, the transfer matrix formulation for a 
specific system can be described and studied. The system studied is a hinged connection 
of two rigid lump masses connected with a linear translational spring and damper 
connected in an infinitely long cable, shown in figure 4-1. The cable is aligned 
horizontally and the tensions are assumed to be constant along the cable. 



R 




Figure 4- 1: Diagram of System 

Where the relationship between the torsional components defined in chapter 3 and the 
linear translational components described here are given by the following; 

Kg = Ka'^ and Rg = Rb^ EQN 4- 1 

An incident wave traveling along the cable can be described by; 

EQN 4- 2 

where A is the amplitude, k the wave number and (o the radial frequency of the wave. 

The subscript / corresponds to the incident wave. 

When the incident wave encounters a discontinuity along the cable, the wave can 
be reflected or transmitted through the region. A comparison of the wave amplitudes 
provides the following; 
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A. 



EQN 4- 3 



7^ = — and = - 
A, A. 



where 7^ is the transmission coefficient and is the reflection coefficient. 

Knowing the rate at which energy is carried by an incident wave along the cable is 
defined by 



1 



2 A 2 



E =—p.Cco' A 
1 ‘ 



EQN 4- 4 



Where is the linear density of the cable, C is the wave speed, ^y is the frequency and A 
is the amplitude of the wave. 

When the wave encounters the system shown in figure 4-1, the energy of the wave 
can be reflected, transmitted, or absorbed into the damper between the two masses. The 
energy of the waves that reflect back from the system can be determined as well as the 
energy of the waves that transmit through the system. The remaining energy, if any, 
represents the energy absorbed by the system, considering that energy must be conserved. 
The ratio of these energies are described by 



and 



therefore 



p - — ^ alternatively p = 

E, 



(±) 



y = — - alternatively y = — 

£, Z, 

£ 

a - — ^ alternatively a ■= \ - p - y 




EQN 4- 5 



EQN 4- 6 



EQN 4- 7 



where the subscript / stands for incident, r stands for reflected, t stands for transmitted, 
and a stands for absorbed. 

But the reflected and transmitted power ratios can also be described in terms of the 
reflection and transmission coefficients (equation 4-3). 
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EQN 4- 8 



^ = |/(,r’ andy = |^|7;|= 
where a is still described as in equation 4-7 . 



4.2 Resonance Frequencies 



The first step of the analysis is to determine the frequency at which the system 
would have the maximum amount of energy absorption. This occurs at a frequency when 
the reflection coefficient is zero, thus all of the energy is either transmitted or absorbed. 
Using the elements of the transfer matrix, the expression for the reflection coefficient is 
provided by: 



R e(o) ) 




Z |Z 2’t i2(“ ) + J “ ■ (z ]'l 22^ “ ) + Z 2 • 1 ]( “ )) 



EQN 4- 9 



Where the elements of the square transfer matrix are tij(a>), for i and j equal to 1 to 2 and 
the impedance of the input and output cable are Z\ and Z 2 , respectively. 

In addition to the reflection coefficient, the transmission coefficient can also be 
expressed in terms of the transfer matrix. The expression for the transmission coefficient is 
provided by; 



T 5 (d) ) =- 



j •2 (i)-Z j 

22( d))-hZ2l 



EQN 4- 10 



Reducing the expression for the reflection coefficient aids in the location its zeros; 



Rc(u) = 



/24'*''22'»2222 S'* \ 

-hd) Z Z 2 “ I Z I* t 22<^W) -( 1 ) Z 2 lii(w) -l-2r(i) Z J i 22 <<o) l 2)(<») ^ 4“2i ( 1 ) Z I Z 2 “ t I ](w) l ]2(<^)/ 

I 2 2 4-*-*22'*22 22 

\l2](0)) ~2l2](d))(i) •ZJZ2l^2fto)-^(l) Zj^Z 2^f j2(^) T(i) Zj^i 22)( ^ ) ■+* 2 O) Z j i 22 ^ (i>)Z2lii(9)J'hd> Z 2 j 



EQN 4- 11 



22 



The numerator of this expression can be solved for values of to such that it equals zero. 
These values of co are the frequencies at which there is no reflection and the system has the 
potential for maximum energy absorption. 

A second method of calculating the frequencies of the zeros of the reflection 
involves the use of the input impedance, defined by; 



z 



m 



and 

velocity 



Zm(0)) 



•wZ2t ii(w) 

2 

-0) -Z 2-t i2<0J ) +} •wt22^“^ 



EQN 4- 12 



By using the principal of impedance matching, the values of co such that the 
imaginary part of this impedance is zero will provide the values of the frequencies where 
the potential for maximum energy absorption will occur. This leads to: 



Imaginary part of z ( o ) = 



-0)^-Z 2 i n(w)-t J2(w)-t2i(«)'0>*t22^“ ^ 

^^2,22,2 
0) 'Z2 ^12^^) ^ 22^ ^ ^ 



EQN 4- 13 

The numerator of this expression can also be solved in order to find the values of to such 
that the expression equals zero. These values are the resonance frequencies of the system. 



Table 4- 1: Input Values for the Determination of the Reflection Coefficient Zeros 



Characteristi 

c 


Value 




Characteristi 

c 


Value 


a 


.0254 m 




b 


.0254 m 


d 


.1524 m 




f 


.1524 m 


c 


.1778 m 




e 


.1778 m 


M, 


4.52 kg 




Mj 


4.52 kg 


Ii 


.035 




I 2 


.035 kgm“ 


Pi 


1.144 kg/m 




P2 


1.144 kg/m 


Tm 


1000 newton 




Ti 


1000 newton 


K 


10 newton/m 




P 


0 newton/(m/sec) 


H 


1000 newton 









An analytic approach and a corresponding numerical solution to this problem is 
outlined in Appendix A. With the value of the damper set to zero, a numerical solution 
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can be obtained for both the zeros of the input impedance expression and the reflection 
coefficient expression. A numerical solution to this problem requires the input of 
characteristic values of the physical parameters of the system. In order to compare the 
simulation with a real world model of a cable system modeled by Vandiver [6], the values 
of the system shown in Figure 4-1 were set in accordance with Table 4-1 . 

The specific values shown in Table 4-1 represent a real world example of some 
oceanographic mooring components. Using these characteristic properties, a computer 
simulation of the system was conducted and the value of the frequencies of the zeros of 
the reflection coefficient were obtained using the secant method. The resulting values 
were Hi = 42.6 rad/sec and Q 2 = 94. 1 rad/sec. Using the impedance method, the same 
values were obtained. The values of the reflection coefficient and the transmission 
coefficient are plotted in Figure 4-2. 




(radians/ see) 

Figure 4- 2: Graph of Reflection and Transmission Coefficients Using Input 

Parameters from Table 4-1 (R=0) 

For the curve shown in figure 4-2, the value of the damper was set at zero to 
observe the location of the minimum values of the reflection coefficient with no energy 
being absorbed in the system to aid in determining a starting value for the secant method 
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analysis. The location of these minimum values is not affected by a change in the damper 
value, R, although the values of the reflection and transmission coefficients are. This is 
shown in Figure 4-3, by increasing the damper to 500 N/(m/sec) The difference 
represents the amount of energy that is absorbed in the system. 




Figure 4- 3: Graph of Reflection and Transmission Coefficients Using Input 
Parameters from Table 4-1 (R=500 N/(m/sec)) 



4.3 Optimum Value for the Damper 

The determination of the optimum value of a damper to absorb the maximum value 
of energy is frequency dependent. The maximum amount of energy will be absorbed at 
the resonance frequencies determined in the previous section. Using the two resonance 
frequencies as a basis, a search was conducted to determine this optimum damper value 
using the same expressions outlined in Appendix B. This time, however, all of the 
equations are dependent upon the value of a particular damper with the frequency held 
constant at either Qi or Q 2 - This investigation emphasizes the higher resonance frequency 
because there is more relative motion between the two masses that leads to more energy 
absorption in the damper as seen in the modal sketch in Appendix A. 
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In order to determine the value of the damper that will lead to the maximum 
amount of energy absorbed for a particular frequency, the absorption power ratio must be 
analyzed. Using equation 4-6, the derivative of the absorption power coefficient with 
respect to R is described by: 

^ - il_iL 

6R 5r 6r EQN 4- 14 



Substituting equation 4-7 into equation 4-13, the following is obtained; 
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EQN 4- 16 



Because the order of the partial derivative operation and the complex conjugate 
calculation are interchangeable, the final expression can be rewritten as: 

«R Ur 'Ur// zi Ur 'Ur// EQN 4- 17 

The partial derivative of the transmission coefficient and the reflection coefficient 
require that the partial derivative of the transfer matrix with respect to R be obtained. 

This analysis is described in Appendix B and is based upon the work of Li [1 ]. The 
analysis provides expressions for the two partial derivatives: 
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and; 



“ ‘(^ 22jnme^*^^'^ 1 ~*~hlj>nine^*^^'^2) (“ 1^ 2 * I2_prirne^ ^ ~ ‘ 21 jOTine^ 

[(*21^^^”*^ '^r^2^12^^^)'*~j 22^ R)'Z I -ht |j(R)Z 2)j 

EQN 4- 19 

The determination of the optimum value of the damper to absorb the maximum 
energy becomes a numerical solution, given the physical characteristics of the mooring 
system. The extreme case is the replacement of the second mass to encounter the 
incoming wave in Figure 4-1 (M 2 ) with a hinge termination. The optimum value of the 
damper would be when the damper impedance was equal to the impedance of the input 
cable: 

L} R Rb ^ 

therefore: Damper Ratio = or Damper Ratio = 7 

b Kpu^un, pcL 

EQN 4- 20 

Where R is the value of the damper, b is the transverse distance from the centerline of the 
system to the location of the damper, L is the length of the first mass ( h in Figure 4-1 ), p 
is the linear density of the cable, and c is the speed of the wave traveling along the cable. 

A comparison with the value of this optimum damper will provide a great deal of insight 
on the performance of the evaluation and design process. 

The design process described in Appendix B was based upon the fact that the user 
will know the value of the linear density of the cable and the range of tension expected in 
the system. From these two parameters, the impedance of the cables (Zi and Z 2 in Figure 
4-1) can be easily determined. Using those impedances as the entry point in the design 
process, a new set of resonance frequencies must be obtained using the process described 
in Appendix A. Based upon the highest of the two frequencies calculated, the optimum 
value of the damper is numerically determined using the secant method for the derivative 
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by substituting equations 4-19 and 4-18 into equation 4-17. This process is done over a 
range of values for the cable impedance, and the results are 




Figure 4- 4: Plot of Optimum Damper Value versus Cable Impedance Using 

Input Parameters from Table 4-1 and Evaluated at the Corresponding ^2 

plotted in Figure 4-4, The range of the cable impedance varied from 10 to 200 kg/sec, 
maintaining the value of the linear density of the cable constant. The linear density as well 
as the rest of the physical characteristics are the same as the values as outlined in Table 4- 
1, with the exception of the tension at Tj and Tj.| which must be altered in order to 
maintain the desired value for the cable impedance maintaining the linear density of the 
cable constant. 

A plot of the values of the absorption power ratio versus the damper values for 
varying increments of the cable impedance is provided in Figure 4-5. When the damper 
values are divided by the cable impedance for the particular variant (Z|=Z 2 =Zc), the 
curves closely follow each other, with a maximum absorption power ratio of 
approximately .5 at a damper ratio of approximately .11. This is shown in Figure 4-6. 

The variance of the cable impedance does not have a drastic affect on the maximum 
absorption power ratio. 

The variation of the optimum damper is also affected by the mass of the 
component masses Mi and M 2 . With the other physical characteristics held constant, the 
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R {ncwu>n/(nv/scc)) 

Zc = 10(kg/sec) 

Zc = 20(kg/sec) 

■ “ Zc = 30(kg/sec) 

“ ’ Zc = 40(kg/sec) 

Figure 4- 5: Plot of Absorption Power Ratio for Varying Input Impedances Using 
Input Parameters from Table 4-1 and Evaluated at the Corresponding Qi 




Damper Ratjo 

Zc = I0(kg/sec) 

Zc = 20( kg/sec) 

■“ Zc = 30{ kg/sec) 

■ ‘ Zc = 40( kg/sec) 

Figure 4- 6: Plot of Absorption Power Ratio Using Input Parameters from Table 4-1 
and Evaluated at the Corresponding n 2 for Varying Input Impedances 
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Damper Ratio 

— MR = 6.48 
MR= 12.96 

— MR = 25.92 
MR = 51.85 



Figure 4- 7: Plot of Absorption Power Ratio for Varying Mass Ratios (Mi=M 2 ) 
Using Input Parameters from Table 4-1 and Evaluated at the Corresponding Qz 



change in the optimum value of the damper with the varying values of the masses(Mi=M 2 ) 
is shown in Figure 4-7. Here the mass of the two components is compared to the mass of 
an equivalent length of cable. The mass ratio is defined by: 



A^ = 



A/, + U, 
2p^d + 2pJ 



EQN 4- 21 



Note that the increase in mass tends to increase the absorption power ratio, and it 
asymptotically approaches a maximum value of .5 when both masses are increased the 
same amount. If the two masses are not altered equally, varying the mass ratio will alter 
the absorption power ratio differently. Figure 4-8 shows that an increase in the mass of 
Ml while maintaining M 2 constant tends to decrease the value of the absorption power 
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— MR =12.96 
MR =19.44 

— MR = 32.41' 

— MR = 58.33 



Damper Ratio 



Figure 4- 8: Plot of Absorption Power Ratio Varying Mass Ratios (Alter Mi only) 
Using Input Parameters from Table 4-1 and Evaluated at the Corresponding CI 2 




— MR =12.96 
MR =19.44 

— MR = 32.41 

— - MR = 58.33 



Damper Ratio 



Figure 4- 9: Plot of Absorption Power Ratio Varying Mass Ratios (Alter M 2 only) 
Using Input Parameters from Table 4-1 and Evaluated at the Corresponding Q 2 
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— Alter Ml Mass Ratio [(M,+M2)/(2p.d+2p201 

Alter M2 



Figure 4- 10: Plot of Absorption Power Ratio varying Mass Ratio (Alter Mi or M 2 only) 
Using Input Parameters from Table 4-1 and Evaluated at the Corresponding Q 2 



ratio, while Figure 4-9 shows that an increase in the mass of M 2 while maintaining Mi 
constant tends to increase the value of the absorption power ratio. Figure 4-10 shows the 
trend towards reaching the ultimate value of total absorption at the second resonance 
frequency as the mass ratio of the system is increased by altering the mass ratio with 
increasing Mi or M 2 only. 

Another physical characteristic of the system involves the lengths of the 
components. If the lengths of the components are changed, the natural frequencies and 
the corresponding values of the optimum damper will also change. Figure 4-1 1 represents 
the effect that changing the length of Mi and M 2 has on the performance of the system, 
using the following relationship for the length ratio; 




EQN 4- 22 



Where d and f are the lengths described in Figure 4- 1 . Note that the location of the center 
peak for the value of the damper in Figure 4-11 is not constant for the five different length 
variations as it was for the mass and the impedance variants. An increase in the length of 
Ml while keeping the length of M 2 constant causes a decrease in the absorption power 
ratio, whereas the increase in length of M 2 while keeping Mi constant increases the 
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Zc = 33.823 kg/sec 



/ 

/ 

/ 

/ / 
I ' 



0 0.5 
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1.5 2 2.5 

Damper Ratio 



Figure 4-11: Plot of Absorption Power Ratio Varying Length Ratio Using Input 
Parameters from Table 4-1 and Evaluated at the Corresponding Hi 



absorption power ratio. The geometry of the system is the dominant term in determining 
the location of this peak value. While altering the length, the other constants outlined in 
Table 4-1 remained constant with the exception of the moment of inertia which is affected 
by the geometrical changes. 

The common factor in both the mass and the length analysis is that an increase in 
the moment of inertia in the first mass to see the incoming wave as defined by the analysis 
(Ml), whether caused by an increase of mass or a geometrical change, decreases the 
absorption power ratio, and less energy is absorbed in the system. An increase in the 
moment of inertia of the second mass in the system to see the wave (M2), will 
correspondingly increase the absorption power ratio and improve the system performance 
if the location of the external force can be identified. 

Even though the lower resonance frequency provides a location when the 
reflection coefficient is zero, the system does not have enough relative motion between the 
two masses to absorb a great deal of energy. This is shown in Figure 4-12. 




Omega 2 = 94. 1 rad/sec 



Figure 4- 12: Comparison of Absorption Power Ratios for the Two Resonance 
Frequencies Using Input Parameters from Table 4-1 

The complexities of the algebra involved with the formulation of the design tools 
could lead to an error that is not readily noticeable. An attempt was made to validate the 
results of this analysis with an independent source. Hayner [7] used a spectral finite 
element method to calculate transmission power, absorption, and reflection in cable 
systems. Similar to conventional finite element methods, a stiffness matrix is formulated, 
which relates the element displacements to the applied forces at nodes due to either other 
elements or external forces. The stiffness matrix includes mass and dissipative effects, as 
well, by letting the matrix elements be complex. Formulating the cable dynamics in terms 
of a stiffness matrix brings to bear well known finite element solution procedures 
discussed by Bathe [8] for solving the response of interconnected elements. The spectral 
finite element differs from a conventional finite element in that the formulation is based on 
a wave solution, which is exact to within the assumptions inherent in the governing 
differential equation. Two additional matrices relating element displacements to wave 
amplitudes and forces to wave amplitudes are used to calculate reflection, transmission, 
and absorption coefficients. The spectral finite element approach may also be considered a 
generalization of the transfer matrix approach discussed by Fertis [9]. 
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Figure 4- 13: Error Comparison with Hayner Method of Transmission Power Ratio 



Using Input Parameters from Table 4-1 and R= 20 kg/sec and Ii=l2=.005 kgm‘ 



A comparison of the two methods is provided in Figure 4-13 and 4-14. The input 
data is provided in Table 4-1 with the exception that the value for the damper is set at 20 




“2 



Figure 4- 14: Error Comparison with Hayner Method of Absorption Power Ratio 
Using Input Parameters from Table 4-1 and R= 20 kg/sec and Ii=l2=.005 kgm^ 
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kg/sec and the moment of inertia of the two masses (I i= I 2 ) is 005 kg m^. Figure 4-13 is 
the absolute value of the difference in the transmission power ratios for the two methods. 
Figure 4-14 is the absolute value of the difference in the absorption power ratio between 
the two methods. Both comparisons reveal that the two methods provide nearly identical 
results for both the absorption power ratio and the transmission power ratio 
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Chapter 5 



Final Results/ Recommendations for Further Research 

The analysis of the moored cable system with the in-line absorber was very 
effective. For varying values of cable impedance, an optimum value of the damper was 
determined. Throughout the entire range of the possible forced vibrations considered 
from vortex induced forces, there was a corresponding value of a damper that could 
mitigate the effect of the motions. This will allow the tuning of the system to help prevent 
the lock-in of the response of the cable with a forcing function near its resonant frequency. 

Throughout the entire effort to determine the value of the damper that would 
provide the maximum amount of power absorption, the magnitude of the absorption 
power ratio was limited to .5 when alterations were made symmetrically with respect to 
the characteristics of the two masses. This matches well with previous attempts to 
maximize the energy absorbed in a cable system by Vandiver and Li [10], although that 
analysis involved the use of a different type of absorber that had three masses instead of 
two and two springs and two dampers. 

When an investigation was performed on the behavior of the system with increased 
cable impedance while maintaining the linear density constant the results indicated a near 
linear relationship between the value of the damper required in order to absorb the 
maximum amount of energy and the tension seen on the cable. As the tension increased, 
the value of the damper required in order to absorb the maximum amount of power also 
increased. A second independent researcher, Hayner, was able to duplicate the results 
obtained in chapter four with a basis of only the input parameters, an indication that the 
process was correct. 

The alterations of the physical characteristics of the masses had a great affect on 
the behavior of the system as a whole. An increase in mass of the system, if divided 
equally between the two component masses, did not provide as much absorption as a 
concentration of the increase in mass in M 2 alone, where an absorption power ratio of 
nearly 1.0 was obtained when the mass of M 2 was 16 times Mi. The limit to that approach 
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is that the differences between the two masses would have other second order affects on 
the system that are not modeled in this analysis. An increase in the mass of Mi alone 
decreased the value of the absorption power ratio, decreasing optimum power absorption. 

The alteration of the component mass lengths also had a significant effect on the 
maximum power absorbed. An increase in the length of Mi decreased the absorption 
power ratio, and an increase in the length of M2 increased the absorption power ratio. 
These trends follow the observations of the mass changes described above. In both 
instances, the changes reflected changes in the moments of inertia of the component 
masses. When the moment of inertia of M2 was increased such that it resembled a fixed 
pin, the power absorption was increased. 

The first area that requires further research is the application of this process to the 
finite case, where the existence of the boundary conditions at the termination of the system 
would alter the results of this infinite cable application. 

The dynamic response of the system to varying currents or changes in speed 
should be included in a further study of the system. 

The inclusion of the bending stiffhess of a cable would provide a more robust 
model with a more detailed wave propagation. 

The construction and testing of a real world absorber would apply the principles of 
this thesis and validate the power absorption methodology. 

Additionally, an investigation should be made into the effects that multiple in-line 
damper systems will have on the behavior of the system as a whole. 
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Appendix A 

Determination of the Zeros of the Reflection Coefficient 
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Figure A-1: Schematic of the System 



Based upon the formulations of the transfer matrix equations from Li, the 
following analysis describes the approach to determining the location of the frequencies 
where the power absorption will be a maximum. A computer program, Mathcad 5.0®, 
was used in the calculation process. 

The first expressions to determine the transfer matrix are expressed in terms of the 
values of the particular dimensions: 
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The Intermediate equations are expressed: 
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The intermediate expressions continue; 
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Now knowing these, the transfer matrix can be determined; 
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Now the reflection and transmission coefficients are determined 
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The Power Reflection, Transmission and Absorption Ratios are; 
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The input impedance is: 
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In order to determine a numerical solution for the analysis, physical characteristics 
of the system must be input. These values used here and throughout the thesis are based 
upon a real world oceanographic mooring described by Vandiver, discussed in chapter 4. 
The following are the physical characteristics of this system: 
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A curve of the Power and Transmission Coefficients provides an indication where to 
search for the Resonance Frequencies of the system 




(raduns/sec) 

Figure A-2: Graph of Reflection and Transmission Coeflicients Using Input 

Parameters from Table 4-1 (R=0) 

Based upon this, it appears as if the best absorption would occur at the frequencies 
at the location of zeros of the Power Reflection Coefficient. Using the secant method, we 
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try to obtain these values at the location of the curves where the Power Reflection 
Coefficient appears to go to zero. 

The reflection coefficient can be reduced to; 



Rc(a)) 



I 24222222222 32 \ 

\^-l2i(0>) -b w Z J Z 2 t 1 ^2^“^ - 0) Z 2 t n(©) -b2r(i) Z j I 22 ^ “ > 21^“^ “ '^1^2 M h2^“^j 

(2 .2 422 222 22 2?2\ 

yl2j(w) — 2t2j((»)>w ■Z|Z2’tj2(f‘J)'b(i) ’Z| Z2't]2!^(i)) +(») "2j t 22 ^ "b 2 ( 1 ) | t 2 ^ (i))Z2lj|((i))+ti) Z-> t|j((t)) j 



And to determine the zeros, we solve for the values of w such that the numerator 
goes to zero using the secant method. 

The first data point is close to 40 rad/sec for the initial guess a =42.564 — 

sec 



The expression for the solution is provided: 
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The resulting resonance frequency is: 



0) : =42.564-— 



This convergence represents an iteration between estimate and calculation until the 
solution converged to the third decimal point. 



The second data point is close to 90 rad/sec for the initial guess 
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rad 



The expression for the solution is provided; 
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The resulting resonance frequency is; u =94.13 — 

see 



This convergence represents an iteration between estimate and calculation until the 
solution converged to the third decimal point. 
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A second determination of the natural frequencies relies on the determination of 
the value of w when the input impedance matches that of the input string. Because the 
string has a purely real impedance, the value of co such that the imaginary part of the input 
impedance is zero will provide this frequency. 

The imaginary part of the input impedance is described by: 



Again, the basis for the determinations of the value of to such that the expression is 
zero is the location of the minimum values of the Reflection coefficient observed from the 
graph, and using the secant method, we can determine the resonance frequencies. 

The first data point is close to 40 rad/sec; u =42,564 — 



This convergence represents an iteration between estimate and calculation until the 
solution converged to the third decimal point. 




see 



The expression for the solution is provided: 




The resulting resonance frequency is: 



0 ) : = 42,564 



The second value appears to be around 90 rad/sec 




The expression for the solution is provided: 




The resulting resonance frequency is: 



(I) = 94.13 



sec 
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This convergence represents an iteration between guess and calculation until the 
solution converged to the third decimal point. For the impedance method, the result was 
not as sensitive to the initial estimate of the value of the frequency, but the Reflection 
Coefficient and the Impedance Method provided the same results. 

A sketch of the two mode shapes that correspond to the location of the reflection 
zeros are provided in Figure A-3. The sketch represents an estimate of how the two mode 
shapes can be visualized, although a solution to the total system motions was not 
performed in order to provide this sketch. Position A and position B are the two extremes 
of the mode shapes. The slight relative rotation between the two masses in mode 1 
represents the amount of absorption that occurs in this mode shape. The extreme relative 
rotation between the two masses in mode 2 represent the maximum amount of absorption 
that can take place in the system at the second resonance frequency. 




Mode 1: Slight Absorption 



Cable Cable 




Mode 2: Maximum Absorption 

Figure A-3: Mode Shape Sketch of System 
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Appendix B 



Determination of the Optimum Damper Values 



Based upon the formulations of the transfer matrix equations from Li, the 
following analysis describes the approach to determining the optimal value of the damper 
operating at the second resonance frequency. A computer program, Mathcad 5.0®, was 
used in the calculation process. The first difference between this formulation and 
Appendix A is the fact that the relationships are dependent upon R, not co. 



Cl |( R) ■ — T j (jiiuyg 'd + H'C-C-K** — (<!)) *I I "i-j ’CiJ'R'b Ol 2( R) *~T — (©) ^ 'Ci) ‘R*b 

a j(R) :=K t^+j -o R b^ 

Intermediate expressions are described; 

-(m j-l + M vf)-a ,(R)-M ,-c-(a j(R))+<i)^-M 

T j(R) ^ 

M 2’®" j’d M 2*^)' 3( R "t" ® *M j*M 2*'?c'd 






T 2<R) •=- 



M2-e-(a ,(R))-I-(m , d+M2'h)-(a jCR)) +( 0^-M ,-M 2 cc-d 



T " 



r 

i 

M J dlc + 


'[(ai(R))-l-h(a3(R))-h]‘ 

2 

(j) 




j^M 2-c-(a j(R)) -h(M j-d-l-M 2‘h)’(a 3<R)) j-M 2-c-rd 



1 4(R) =- 



l*l 4“ M ‘ 2(R)j'bMj*C‘Ct2(R)"t"® *Mj'M2‘®’c*f 



(tt2(R)h + 0 3(R)l) / , 

^ L - IM , d-r + M 2-f^ h 
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The intermediate expressions continue; 

(a 2 (R)h + 0 3 (R)l)’ 






+ M 2 cfh 






2 2 

-a-JR) (M , d-f-M 2 h)-M 2 e-(l 3 (R)+ 0 ) M . M yc d 
T 6(R) •= 7 X 



Now knowing these, the transfer matrix can be determined; 



T 5(R)+t 2(R)T 6(R) 

1-7 3(R)-7 6<R) 

_T 4<R) + 7 ,(R)T 6(R) 



t ,2(R) =T 2(R) + 7 3(R) 
t 2 l(R) 



tii(R) =7 i(R)-t-7 3(R)- 



7 4(R)-t-7 i(R)-7 6<R) 



'-7 3(R)-7 6<R) 



1 s(R)-I -7 2 <R )-7 6 (R) 
t 22 (R) — 



1-7 3(R)'7 6<R) 1-7 3(R)-7 g(R) 

Now the reflection coefficient is determined 
^ <»-(Z i-t22(R)-Z2‘ n(R)) 

l2l(R)-® 'Z j'Z 2't |2<R)+j '“(Z jt 22(R) +Z 2't j j(R)) 

The transmission coefficient is determined: 



T,,(R) 



j J 



t ~ ® 2 ^ 12 ^^) I I 22 ^ R) + ^ 2 '^ 

The impedance is: 



Zm(R) = 



t2i(R)7-j (i) Z2t ii(R) 



"<D *Z 2 ‘lj 2 (R)'t~j ’W’t 22 ^^^ 



The Power Reflection, Transmission and Absorption Ratios are calculated 

a(R) ; = i-P(R)-7(R) 



P(R) ;=R^(R) R^(R) 



7(R) T.(R) T.(R) 
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Now the determination of the optimum value of the damper value to obtain the most 
absorption can be done by taking the derivative of the Power Absorption Ratio with 
respect to R. 



^ 

5R 5R 5R 



Because 



a = 1-P-Tf 



knowing 



=R(.R 



leads to: 



5R 



_-8(rcrJ 



5r 



5R 



or 



l&R 5rJ 


! 1 


^5T . — 5T J 


R -t-R ^ 

\ 5R 6R / 


' z\ ' 


i c 

Ur 5r / 



-^2 — 
T- T -T 



5a _ 
5 r 



This derivative requires the calculation of the derivative of the transfer matrix from the 
constituent components of the system. The intermediate equations are provided. 

T 1 i*(l^c)-hM2'f]‘[M2-e-(a i(R))+ (m j-d-HM 2'h)-(a 3(R)) j-M2*c*<rd] 

lj>riinc_b<*^> i*d-hM2-(c+h)]-[-(M j-l-HM2‘f)-a j( R) - M j-c- (a 3( R)) j-M 2'f c^] 






j 1 _prime_.(R)+T l_prin,e_b( R>) 



2_primc 



M2*e-^a j(R))-h(^M j-d +-M 2*h)- (a 3 < R)) +w -M j-M 2'C-C'd 
^(R) •=-|-^j*[M2-e-(a j(R)j-t*(M j-d + M 2'h)- (a 3( R)) -hco^-M j-M2*c-e*dj 



^ 2_primc_b<^> I-d + M2-(e-t-h)]' 



2_primc 



(R): 



j -(O-b '(y 2_pnmc 2_primc_b^^^) 

I^M 2’®" 1^ 1*^ d“M 2'^)’ 3^ R)) 2’**’®'^ 
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• (a j(R)) -I- (m j-d-i-M 2'hy (^a 3(R)) j-M 



3 j,rimc_b<R)- = [M i-d4-M2-(c^h)]. 



M j d I c-h 



[(a i(R))-l-h(a3(R))-h] 

2 

b) 



3_primc 



(R):=- 



-j ’Wb 3_prime_b^ 

[[-M 2 'e-(a |(R)) - (m ,-d + M j-h)- (a 3 (R))] - a^-M ,-M j'ce-d 



^ 4_primc_a^^^ 



'[m j-(l4-c)-fM2*f]* 



(a2(R)-h-f-a 3(R)-|) / ~ ^ \ 

(^M j-d-r-hM2'r-h] 



"f 4_priroc_b^^^ 



hjH 

2 I 

© / 



M 



J-I 4-M 3(^)) 2’®*'^^ 



4_priinc 



(R) : = 



. _ 'j ® ^ 4_primc_a^^^ 4_priinc_b^^^) 



(°^.'‘,''"“».'‘,")-(m,.ii^m,i^i,) 



L-i 

J •W'b 



5_priinc 



(R)-- 



h-H 



(a 2(R)‘b -ha 3(R)1) / 2 

; 



d I N-M2' 



f.)- 



2( R)’h -ha ^{R)-!^ 



11 

+ M2-c-fhi 



(a2(R) h-ha3(R)-l) / ^ ^ \ 

(^M i-d-r-hM 2-r-hj 



^ 6_priinc_a^^^ “ 



- j*d -|-M 2’(c-hh)j’ 



(a 2(R)-h-ha 3(R)1) / ^ 2 \ 

[M j-d-l -hM2-r-hj 



^ 6_prunc_b^^^ 




-a 2( R)' j'd “hM — M 2’^ct 3( R) "hb 



•M j M 



2 . 
2‘e d 



•- ^ ^ ^ 6_primc_a^ 6_pfimc_b^ 



(a2(R)-h-+-a3(R)*l) 



(m ,-d- 



d I N-Mj-rh 
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The elements of the derivative of the transfer matrix are then calculated; 



^ 21_pnmc _.(R) =T 4 _primeCR> + T l_pnn.e(R)'T 6<R)n 1<R)T 6 _pnn«(R) + T 3 _jmme(R>T 4<R)T 6<R) 

*21 _prm«yR)'^-T 4 6 <R) + T 4 <R)T 6 

*(t i(R)T 3 _prim*(R)-T 1 j,rime(R)l' 3(R)) (t 6<R>)^ 



„ ,_‘21j)riine a(’^)^"‘21_priine_b(’^^ 

‘ 21_prime(R) ^ -T 

(i-r 3 (R)t 6(R)) 



^ 22 _priine_a^^^ ■ t Sjxime^^^ ‘*' 1 ' 2_pfiine^*^^ t 6 {R)+T 6_priine(R) "•'T 3 _prime^ 5 (R) T fi(R) 



‘ 22 _pmneyR)-=-T 5 jpnmc^^^l 6 <R)+T 3 <R)T 5 (R)T 6 _prime(R)- 

^ (r 2(R)1' 3 j,rin,e(R)-T 2 j,rin.e(R)T 3(R)) (t 6<R))^ 



^ 22_primc^^^ 



^ 22_primc_a^ R) -t-t 22_priine_b^'^^ 

(l-T 3(R) T 6(R))^ 



‘ ll_ptime(^) "T l_prime(’^)+t 3_pnme('^> ‘21< 3<'^> ‘ 21_priine('^> 

‘ 12 _priine(’^> "t 2_prime(’^) + 1' 3_prime( ‘ 22('^) ‘‘‘I' 22_priioe(^> 

The calculation of the derivatives of the Reflection and Transmission Coefficients with 
respect to R are provided: 



dRc_dR_j(R) ^-0) 



dRc_dR_b(R) ;-j 



2 12 j)rintc(^^’^ 21^ 1 l_prinic(^^ ^ 22 ^^^ ' ^ 22 _primc (^0 

i2_primc(R)) " 21 22( ‘ 2l(R> ‘ 22_prune(R)) 



dRc_dR(R) :=- 



2Z j (dRc_dR_a(R)ri-dRc_dR_b(R)) 



2j( R) — 2 ^ 12^ 22^ R)’Z j -h t 1 j( R)'^ 2 ) ] 



$Rc 

6R 



dTc_dR(R) =2Z 



2 / 

® * (5 22 prime' 


(R)'Zl + ‘lI prime(R)Z2)+j ■' 


V" 1^ 2 ‘ 12 primeC^^ ® ‘ 


21 prime^*^^) 




’1 


i^t2i(R)“W ’Z j-Z 2't j2^ R)) j 


22 ^^^'^ 1 2 )] 


r 



6Tc 

5R 



The operation of the complex conjugate and the partial derivatives derivative are 
interchangeable, which leads to: 

6R 



|^dRc_dR( R ) R j,( R ) + R R )• dRc_dR( R )) (dTc_dR( R )• T J R ) ^ T R) dTc_dR( R )) 
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Again, in order to complete a numerical analysis, a set of physical characteristics of 
the system must be input. In the case of this example, the previous values from Table 4-1 
are used again with a concentration at the second resonance frequency of the baseline 
system. Evaluated at the higher of the two resonance frequencies: 



a =.0254 m 


d ;=.1524 m 


c.=. 1778 m 


h -c-f-d 


M J :=4.52 k 


b :=.0254 m 


f :=.1524 m 


c =.1778 m 


1 =c-|- 


M 2 —4.52 k 


P , ;=i.i44.1« 
m 


p,. = 1.144-1* 


kg 

7 —’ll #771. S, 


7 * —'I’l 1- 


newton 

V in 


sec 






m 


ICC 


m 



_ 1 _ 2 

H: = 1000-newtoo T j Tj: 

P 1 P2 



1 ,.= 



M , (2d)" 



12 






' 2 - 



12 



. ^ 

w .=94.13 

sec 



In order to obtain an initial estimate of the value of the damper for maximum 
absorption, a graph of the Power Adsorption Ratio is provided. This will provide the 
range for the initial estimate 




Figure B - 1: Plot of Absorption Power Ratio For Varying Values of R Using The 

Input Characteristics of Table 4-1 
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The expression to determine the optimum value for the damper is provided by: 
The initial estimate is around 500 newton/(m/sec) from the graph above 



R : 



=464.343 



nev^ton 



m 

see 



Bc«t_R ; = 


root 


- (dRc_dR(R) R g(R) + R ^(R) dRc_dR< R)j ... 






+ - 


Z 2 / 

— (^dTcJRCRVT e(R)+T j(R) dTc_dR(R)j 






__ newton 

Bc«t_R =464.343--; — p 

fe) 

This represents an iteration until the result was within three decimal places of the estimate. 
An iteration was performed until the solution was within three decimal places of 

the initial estimate. The value of the optimum damper for this case is: 464.343 kg/sec. 
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